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I. INTRODUCTION 

Quantum many-body systems often can be simplified 
considerably by choosing an appropriate basis. This is 
the first step where the qualitative understanding of the 
system enters. If the problem is first formulated in a ba- 
sis which is not optimum for the subsequent treatment 
it is reasonable to change the basis. This is a unitary 
transformation. Famous examples of such changes of ba- 
sis are the fermionic one that leads to the wave function 
of the superconducting BCS phase or the bosonic Bogoli- 
ubov transformation which leads to the spin waves in a 
Heisenberg quantum antiferromagnet. These examples 
refer to mean-field treatments. The complete solution of 
a many-body problem by a simple analytical transforma- 
tion is rarely possible. 

Ten years ago, the proposal was made to use unitary 
transformations beyond the single particle leveli<2iii. The 
unitary transformation is performed in a continuous fash- 
ion whence the name continuous unitary transformation 
(CUT). This means that it is parametrized by a contin- 
uous parameter I running from to oo. At I = the 
Hamiltonian is unchanged and at Z = oo it is in its final 
form. The fundamental idea is that the transformation 
at a later stage, i.e., at larger values of I, can take into ac- 
count the changes introduced already at an earlier stage. 
The transformation can be set up such that processes at 
larger energy are treated before those at lower energies. 
Then the transformation has properties similar to Wil- 
son's renormalization group approach^. The objective is 
of course to simplify the model by the transformation. 
At best, the final Hamiltonian has become diagonal. 

The aim of the present paper is to illustrate the CUT 
approach, and to study it for a simple model where exact 
(numerical) data is also easily accessible. The simplicity 
of the model permits to present the equations on the basic 
level and also to go to higher orders. Our presentation 
is intended to serve as introduction to the CUT method. 
This will enable the interested reader to generalize the 
approach to more elaborate models. 

The model of our choice is the quantum anharmonic 
oscillator which is defined by the usual quadratic ki- 



netic energy and a potential with quadratic and quartic 
parts. We will abbreviate the oscillator comprising both 
a quadratic and a quartic part in the potential by QO; 
the one which has only a quartic one is called the pure 
quartic oscillator (PQO). This model has often been in- 
vestigated as testbed for more complicated situation, e.g., 
it was used to test other renormalization treatments^. 

The bottom line of our investigation is that the CUT 
method works extremely well providing very good quanti- 
tative results. The starting point of the CUTs can easily 
be optimized by incorporating prior understanding of the 
physics of the underlying model. Then only very moder- 
ate effort suffices to obtain satisfactory results. 

Following this Introduction, we present basic facts and 
insights about the CUTs. In Sect. Ill results for the 
energetically lowest states of the QO and the PQO will 
be shown for various ways to conduct the transformation. 
Spectral weights will also be computed in this section. 
In Sect. IV we discuss in which way the CUT approach 
can be systematically improved. Sect. V concludes the 
article. 



II. BASICS ABOUT THE CUTS 

Before addressing the anharmonic oscillator, it is 
worthwhile to study the harmonic oscillator in order to 
become familiar with the CUTs. 



A. Method 

A given Hamiltonian H can be diagonalized by a suit- 
able unitary transform U (UU^ = I). But it is usually 
a very hard task to find this transform U. The idea of 
Wegneri (and independently of Glazek and Wilson 2 * 3 .) 
was to diagonalize the Hamiltonian in a continuous way 
starting from the original Hamiltonian H = H(0). A 
flowing Hamiltonian is defined by H(l) = U'(l)HU(l) 
depending on the parameter I such that H(l — oo) is 
decisively simpler. 

The continuous transform is equivalent to performing 
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infinitely many infinitesimal unitary transforms e~ v ^ dl 
with the generator rj(l) = —U^{l)diU(l). The flow equa- 
tion 



d l H(l) = [ V (l),H(l)} 



(1) 



defines the change of the Hamiltonian. Of course, the 
whole transformation U(oo) = £exp(— J r](l)dl) can 
also be given in one formal expression. Note, however, 
that this requires the use of the Z-ordering operator C 
which sorts from left to right according to ascending val- 
ues of I. This means that the whole transformation op- 
erator can be extremely complicated. For simple models, 
however, the explicit transformation can also be written 
down. 

The crucial point is to choose the generator r\ which 
leads to a simplification of the Hamiltonian. Wegner pro- 
posed to take the commutator between the diagonal part 
.ffd and the non-diagonal part i? n d in a given basis. This 
means that the generator reads r/ w (l) — [H d (l), H n d(l)] = 
[Ha(l), H(l)}. In the basis of the eigen states of H d the 
generator thus reads 

r Htj {l)=(e i {l)-e j {l))H id {l) . (2) 

For finite matrices, this choice was proveni to achieve 

[H A {oo),H(po)] = . (3) 

If -ffd(co) is non-degenerate this equation implies that 
H(oo) is diagonal. Otherwise it implies block-diagonality 
with respects to the degenerate subspaces of -ffd(oo). An 
extension of this proof to infinite matrices is given in 
Appendix ^ Note that there is an enormous freedom in 
what one considers to be diagonal and non-diagonal. For 
practical calculations the only restriction on H d should 
be that it can easily be diagonalized. 

For Hamiltonians of the form of band matrices, Mielke 
proposed another generator involving the sign of the dif- 
ference of the vector indices (see below). This generator 
preserves the band structure during the flow 7 which is not 
the case for Wegner's generator. Independently, Knetter 
and Uhrig observed that a generator based on the sign of 
the difference of the particle number allows to construct 
particle-number conserving effective Hamiltonians&&. If 
Q is the operator counting the number of elementary ex- 
citations and the matrix elements of rj in the eigen basis 
of Q are chosen to be 



(4) 



and the final Hamiltonian satisfies [Q, H(oo)] — 0. Fur- 
thermore, one can show that this type of generator sorts 
the eigen values in ascending order of the particle number 
of the corresponding eigen vectors 7 ^. Analogous gener- 
ators were also used by StemiS*ii for models where the 
sign-function is not required since only one finite absolute 
value of Aq := qi — qj occurs. For instance, if Aq = 0, ±2 



the sign- function only changes the scale of I. But for 
Aq = 0, ±2, ±4 the sign is essential. 

In the sequel, the generators of Wegner's type will be 
denoted by the subscript W; the generators characterized 
by the sign as introduced by Mielke, Knetter and Uhrig 
will be denoted by the subscript MKU. Any coefficient g 
without argument appearing in the Hamiltonian or in the 
generator stands for the running coefficient g(l). Their 
initial values are denoted by (bare) and their final 
values by g^ (renormalized) . 



B. Two Solvable Models 

We illustrate the above for easily solvable models, 
namely a two-level system and a harmonic oscillator. The 
Hamiltonians read 



H 2 = EI-~a z +~a*, 
Hho = E1 + wala + ^ (a) 2 + a 2 



(5a) 
(5b) 



where a 1 are Pauli matrices, and a and are bosonic 
annihilation and creation operators satisfying the com- 
mutation relation [a, cr] = 1. 




00 









FIG. 1: Flows for the two-level system. Panel (a) depicts 
the flow from the MKU-generator; Panel (b) the one from 
Wegner's generator. Thick lines represent lines of fixed points; 
the dots represent the fixed points of the particular circular 
trajectory shown. The arrows indicate the direction of the 
flows for increasing /. 

Let us first consider the two-level system. The count- 
ing operator is Q = (I — c 2 )/2. Hence, the MKU- 
generator is given by !) M8U = —{e/2)ia y . The flow equa- 
tions read 



diE 
dioj 
die 







-cue . 



(6a) 
(6b) 
(6c) 



The diagonal part of the Hamiltonian is the one pro- 
portional to ui. Thus Wegner's generator is found to 
equal i] w = [H d ,H nd ] = [-(ui/2)a z , (e/2)a x ] = ljt) mku . 
So the flow equations are the same as above if we rescale 
the continuous parameter I — > I' with dl' — u>(l)dl and 



3 



= 0) = 0. The flow trajectories in (oj,e) space are 
circles since oj 2 + e 2 = R 2 is a constant of the flow, see 
Fig. H Thus the diagonalized Hamiltonian has the coef- 

7 and £( R ) = E( B K 



ficients oj^ = y/R = y/tJFr 



(By 



FIG. 2: Flows for the bosonic oscillator, (a) MKU-generator 
and (b) Wegner's generator. Thick lines represent lines of 
fixed points, and the dots represent the fixed point reached 
for the particular hyperbolic trajectory shown. The arrows 
indicate the direction of the flows for growing I. 

Second, we consider the bosonic oscillator for which the 
counting operator is Q = a' a. So the MKU-generator is 
??mku = d/2{a* — a ) and the flow equations read 

diE = -d 2 (7a) 
diuj = -2d 2 (7b) 
did = -2ud . (7c) 

The canonical choice of the diagonal part is i?<j = oja^a. 
Thus we have again rj w = [uxva, d/2{a' + a )] = ojr] MKU 
so that the flow trajectories of both generators are the 
same. The expression lu 2 — d 2 is a constant of the flow so 
that the trajectories are hyperbolae as shown in Fig. [21 
For a hermitian Hamiltonian bounded below, i.e., oj > 
\d\ (see below), the final coefficients of the diagonalized 

Hamiltonian read oj ( R ) = Voj( b ) 2 -d(B) 2 and = 
E^ + (oj^ b ' > —oj^)/2. The energy is found by integration 
of d t E = dioj/2. 

As expected from the rescaling dl — ► d^ = ojdl the 
topologies of both flows are the same. The directions 
are the same for positive values of oj. Only for negative 
values of oj the directions differ. The coincidence of the 
topologies of the flows for both choices of the generator 
is special to the simple systems under study here. It 
is not a general feature for more complicated, generic 
Hamiltonians. 



C. Remarks 

First, we want to stress that the CUTs exactly diag- 
onalize these simple Hamiltonians. Note that for the 
harmonic oscillator, the CUTs are infinitesimal Bogoli- 
ubov transformations whose total effect corresponds to 
the single-step Bogoliubov transformation diagonalizing 



the Hamiltonian. 12 This was to be expected since there 
is only one unique way, except for phases, to diagonalizc 
non-degenerate problems. 

Second, we see that Wegner's generator always leads to 
a fixed point. This does not need to correspond to a final 
diagonal Hamiltonian when there are degeneracies as it 
is the case for oj = 0. In contrast the MKU-generator 
is not sensitive to degeneracies, but it does not always 
lead to a fixed point. A fixed direction may occur as 
seen in Fig. |2Ia)i22i The flow in Fig. |2Ja) is similar to 
the Kosterlitz-Thouless flow. The fixed directions are 
reached when the Hamiltonian is unbounded below, i.e. 
oj < \d\. In this case the proof of convergence indeed 
failsiSi A possible cure to this problem is to diagonalize 
—H, or what amounts to the same, to revert the flows. 12 
Physically, Hamiltonians unbounded below occur only for 
open physical systems whose treatment requires special 
care. This is beyond the scope of the present investi- 
gation. It is a general feature that Wegner's generator 
always leads to a fixed point and that it is stopped by 
degeneracies. In contrast the MKU generator is insensi- 
tive to degeneracies, but it may fail to converge to a fixed 
point. 



III. THE QUARTIC OSCILLATOR 
A. Model and Notation 

The Hamiltonian of the anharmonic oscillator under 
study is in first quantization 

f 
2 

where A is a positive constant. We will consider two 
cases. For oj — 0, we consider the pure quartic oscillator 
(PQO) for which a rescaling X — > X/a, P — > aP and 
H -> H/a 2 leads to a unique Hamiltonian. Choosing 
a 6 = 6A it reads 



H 



(P 2 +oj 2 X 2 ) + XX 4 



(8) 



H 



1 



-P 2 



-X 4 



(9) 



For oj ^ 0, we use another rescaling to bring the Hamil- 
tonian in the form 



H 



p2 



X' 



XX 4 



(10) 



The Hamiltonian (|10|l will be called QO in the sequel. 

Both Hamiltonians have been studied numerically, see 
for instance Ref.yji Nowadays, they can be diagonalized 
numerically on a workstation in seconds, so that we can 
easily check the quality of our approximations. 

For both Hamiltonians, we introduce annihilation and 
creation operators 



a(a) 



aUa) = -4= I — -io-P 
\/2 W 



(11a) 
(lib) 
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where a is a free parameter which will be determined 
later. Recall that the QO reduces to the conventional 
harmonic oscillator for A = 0. It is diagonalized by us- 
ing (7 = 1 in Eq. 111|) . The Hamiltonians will be used 
in their normal-ordered form. Normal-ordering is done 
with respect to the bosonic vacuum which depends on 
the value of a. The operators are put to the left of the 
a operators. Thus the general form of the Hamiltonians 
reads 



H = .90,0 + .go,i flta + go^a^a 2 (12) 
+a f2 (52,0 + 52,ia f a) + (52,0 + 92,1^ a) a 2 

t 4 4 

+a 34,o + 54,oa ■ 



For brevity of the notation, we do not denote the explicit 
dependence of the couplings and operators on a in the 
above equation. The coupling constants gty have two 
indices which indicate the powers of o' and a in the cor- 
responding operator. For instance, g^j is associated to 

the operator ■ a° and to its hermitian conjugate. 

For illustration we consider the harmonic oscillator 
(|5b(l . The ground state energy is given by E = <?o,o, 
the single-particle energy by u> = go,ii an d the Bogoli- 
ubov term by d/2 = 52,0- For the PQO, the bare values 
of the coupling constants are 



(13a) 
(13b) 
(13c) 
(13d) 

(13e) 
(13f) 
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For the QO, they are found to be 



7 f \+a 2 + 3Xa 

4 V (T 2 



0SV) 

<?foV) 
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cr 2 + 6A<r 



a + 6Aer 



(14a) 
(14b) 
(14c) 
(14d) 

(14e) 
(14f) 



B. Perturbative Results 

We start by analyzing the QO by perturbative 
schemes. The results shall serve as reference for the fol- 
lowing more sophisticated approaches. 

The first approach is to expand perturbatively in A 
around the diagonal Hamiltonian for a — I. It is 
known that the resulting series is not convergent but 
asymptotic^. Yet the results may approximate the ex- 
act ones if they are restricted to low orders in A. In this 
sense, we computed the first three renormalized eigen en- 



ergies Eq , Ei and E^' to second order in A 



(R) 



bare values are E n = 



E. 



(B) 



(B) 
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■9o,o ■ 



E[ B) = 



o (B) 
9 0.0 
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9o,i 
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and 



2^. 



From the approximate first 



three eigen energies we deduce the renormalized diago- 
nal coefficients of the Hamiltonian g^Q , gjy^ , g^ ■ Fig- El 
displays the relative errors of these coefficients compared 



(R) 

9o,i 



(ex) 

9o,i 



(ex) 



The 



to the exact values, i.e., r% 

curves are labelled 20-PT for second order perturbation 
theory. As expected the errors grow very quickly on in- 
creasing A. 

Instead of passing tediously to higher orders of A we 
choose another route. We consider all terms which pre- 
serve the number of quanta as diagonal. This means all 
terms with equal number of creation and annihilation op- 
erators are kept in the diagonal Hamiltonian. Then we 
perform second order perturbation theory around this 
diagonal part in all the off-diagonal contributions, i.e., 
the terms which change the number of bosonic quanta. 
For concreteness, we give the resulting expression for the 
ground-state energy 



( R )^ <l 

00,0 M 



(15) 



4(sEV)+3< 2 V)) 

Now, we exploit the freedom to choose a. Five choices 
are studied: 

• cthf = 1. This choice diagonalizes the harmonic 
part of the Hamiltonian l|lUfl • Hence we denote this 
scheme 20-HF since it represents second order per- 
turbation around the Hartree-Fock Hamiltonian. 

• a is chosen to fulfill 9<r5oo (°0 — 0, i.e., a station- 
ary point of the bare approximate ground state en- 
ergy is chosen. It can be checked that this is the 
same as to have 32.0 (°0 = vanish. So the bilin- 
ear part of the whole bare Hamiltonian is diago- 

(B) 

nal. This scheme is denoted by 20-VAR-<?q since 
it represents second order perturbation around the 

(B) 

stationary point of g$ q . Conventionally, this choice 
is called second order perturbation around the self- 
consistent Hartree-Fock Hamiltonian. 



(B), 
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FIG. 3: Relative errors ro, ri and r 2 of the first three diagonal 
coefficients g^g, g^i , g^2 (see text) for five different com- 
putation schemes. For the second order perturbation theory 
(20-PT) n reaches a maximum value ~ 23% for A ~ 0.5 
(outside of the graph) and then decreases again; r 2 reaches a 
maximum values ~ 300% for A ~ f.f. For the second order 
calculation around the Hartree-Fock Hamiltonian (20-HF), 
all Ti grow very fast. The results of the second order cal- 
culation around the optimum bare value of the ground state 
energy (20-VAR-(?g 1 ^) and the one around the optimum value 

of the renormalized gj^g (20-VAR-g^Q) are always very close 
to each other. For ro they cannot be distinguished. 



• In analogy to the previous choice we require 
that the bare single-particle energy is stationary 

d i9o B i = °- This choice is called 20-VAR-fl^. 

• Instead of looking for stationary points of the bare 
parameters, we can look for the stationary points of 
renormalized parameters. We denote by 20-VAR- 
(7g q the second order perturbation around the sta- 
tionary point of the renormalized ground state en- 



ergy d a g ( Q J(a) = 0. 

• We call 20-VAR-g^ the analogous variation of the 
renormalized single-particle energy. 

On the mean-field level it is well-established to look 
for stationarity in order to define the optimum choice 
of parameters. To do the same for higher order results, 
here second order, is justified by the argument that in an 
exact calculation the result does not depend on a. The 
value of a only defines a starting point of the calcula- 
tion. If the calculation is exact it does not depend on the 
starting point. Hence it is stationary with respect to the 
variation of the parameters defining the starting point. 
An approximate treatment, however, will not show this 
independence on the starting point. The best that can 
be achieved is stationaritjii^ which therefore defines the 
optimum starting point. 

The results obtained for these five choices are repre- 
sented in Fig. [3J Clearly, the 20-HF result leads to 
no improvement compared to the 20-PT result. The 
four other choices using a variational criterion lead to 

much better results. Apart from the ground state energy 
(ft) 

where the 20-VAR-gQ scheme is the best, the overall 

best choice is the 20-VAR-t/Q { scheme, especially for r%. 
Note that there is no big difference between 20-VAR- 
g$ and 20-VAR-g^. On average, the 20-VAR-^ 

scheme computes the g^f coefficients with an accuracy 
ranging from half a percent to a few percents. Thus we 
conclude from these calculations that the use of a vari- 
ational criterion is always very rewarding. We will also 
use such a criterion to choose the starting point of the 
CUT approaches. 

Finally, we address the PQO without any harmonic 
part. It represents the A — > oo limit of the QO con- 
sidered before. Of course, neither a direct perturbative 
approach nor a Hartree-Fock approach is defined. But 
all the other approaches can be used for the PQO with- 
out conceptual difficulty. The results are given in Table 
[I] Again, it turns out that the stationary point of the 





10 J r 


10*7-1 


10"r 2 


20-VAR-ff^ 


1.94 


4.47 


54.5 


2O-VAR- ft ( ) H ) 


1.94 


4.57 


55.3 


2O-VAR-g B 1 ) 


4.01 


0.59 


-4.07 


20-VAR-^ 


3.45 


0.47 


3.08 



TABLE I: Comparison of the four variational second order 
calculations for the purely quartic oscillator. 

single-particle energy leads to the best results. The er- 
ror in the interaction value analyzed by r^ is reduced by 
a factor of 10 compared to the stationary point of the 
ground state energy. The deviation in the ground state 
energy is smallest at its stationary points, but only dou- 
bled when the single-particle energy is made stationary. 
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C. Non-Perturbative Results 

1. Spectrum 

Let us consider the Hamiltonian given in Eq. (|12fl . Per- 
forming an infinitesimal unitary transformation gener- 
ates some new terms in the Hamiltonian which were not 
present before. Focusing, however, on the states at low 
energies the new terms with at least six bosonic opera- 
tors will be of minor importance. They represent pro- 
cesses which involve at least three particles. Hence we 
discard the new terms for the moment once they have 
been brought into normal-ordered form. Recall that on 
normal-ordering higher order terms, terms of lower order 
occur. These are kept since they matter for the low-lying 
states. 

If the flowing Hamiltonian is truncated in this way it 
stays in the form of Eq. I|12l) . It is natural to view all 
terms conserving the number of bosons, namely 30,0 + 

go.ia'a + go^a^a, 2 , as the diagonal part Hd- The non- 
conserving term make up the non-diagonal part of the 
Hamiltonian. With this choice, Wegner's generator ?y w = 
[Hd,H n d] is computed. It is found to comprise terms of 

the form a 2 — a) 2 ^ and a^a — a^a 5 which are not 
present in the initial Hamiltonian. Hence we omit them 
also in the generator— which reads 



V = V2.o(a -a )+7?2,i( a a — a'a ) (16) 

4.4 4 
+V4,o( a - a )> 

where for Wegner's generator 

V2,o = 2 (g Q ,i + 50,2)32,0 (17a) 
V2,i = 2(50,1 + 3.90,2)32,1+430,232,0 (17b) 
774,0 = 4(30,1 + 330,2)34,0 • (17c) 

For the MKU-generator the coefficients can be read off 
from the Hamiltonian 

?72,o = 32,o (18a) 
V2,i = 32,1 (18b) 
»74,o = 34,o • (18c) 

It is one of the advantages of the MKU-generator that it 
allways has a simpler form than Wegner's generator, as 
exemplified here for the quartic oscillator. This was not 
obvious for the quadratic oscillator or the two-level sys- 
tem which are too simple to exhibit this general feature. 



The flow equations read 

<9;3o,o = -4?72,o32,o ~ 48?74,o34,o (19a) 
<9z3o,i = -4?72,o(232,o + 33 2 ,i) - 12*72,1(32,0 + 32,1) 

-1927,4,034,0 (19b) 
<9;3o,2 = -12772,0.92,1-12772,1(32,0 + 332,1) (19c) 

- 144/74,034,0 
<9;32,o = -2772,0(30,1 + 3o,2 + 634,0) - 24r7 2 ,i34,o 

-12774,0(32,0 + 232,1) (19d) 
<9;32,i = -4772,0(30.2 + 2.94, 0) - 4774,0(252,0 + 932,1) 

-2772,1(30,1 + 330,2 + 1834,0) (19e) 
<9;34,o = -2772.032,1 + 2772,132,0 (19f) 

-4774,0(30,1 + 330,2) • 

The derivation of these expressions is presented in Ap- 
pendix IbI 

We start with the results for the PQO. All the results 
for the eight possible computational schemes are collected 
in TablclTTl The eight schemes are made up from two pos- 
sible generators (W or MKU) and from the four choices 
which quantity is made stationary, see previous section. 

We have to mention that 3,3 and g$ -[ display both a 
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W-VAR- 5 ^ 


0.127 


0.000876 


-15.8 


W-VAR-<C 


-0.0815 


-0.195 


-1.86 


W-VAR- ff ^ 


-0.0612 


-0.279 


-4.83 


W-VAR-3^ 


-0.0426 


-0.287 


-6.22 


MKU- VAR- 


0.0492 


-0.136 


-16.5 


MKU-VAR- g( ( ) " ) 


-0.0608 


-0.232 


-0.951 


MKU-VAR-5,™ 


-0.0555 


-0.289 


-3.51 


MKU-VAR-5™ 


-0.0352 


-0.315 


-6.97 



TABLE II: Comparison of the four variational CUT calcula- 
tions for the PQO with Wegner's or MKU generator. 

local minimum and a local maximum, i.e., there are two 
choices for the stationary point. Here only the local min- 
ima are considered. 

The comparison of Tables [I] and [H] reveals that the 
CUT results are always much more accurate than the 
second order calculations performed in Sec. IIII Bl This 
is of course expected and does not depend on the chosen 
variational scheme. For instance, the ground state en- 
ergy is mostly given to an accuracy of less than 0.1%, to 
be compared to the 2% of the second order calculations. 
Yet the calculation is straightforward. Only six differen- 
tial equations need to be solved which can be done by 
any computer algrebra program. This evidences that the 
CUTs lead to a significant improvement. 

To vary the renormalized quantities instead of the bare 
ones modifies the CUT results more than it modified the 
second order results in Table [I] In the CUT calculations, 
the stationary points of the ground state energy provide 
the best results. This is in contrast to the finding in 
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the perturbative results where the stationary points of 
the single-particle energy were optimum. We stress that 
the variation of the renormalized quantities is superior 
to the variation of the bares one if the focus is set not 
only on the ground state energy but in particular to the 
single-particle and two-particle energies. 

In comparison, the results for the two choices of the 
generator are very similar. The MKU-generator works a 
little better, for instance for the variation of the renor- 
malized ground state energy. We attribute this to the 
fact that it avoids the appearance of terms which create 
or annihilate more than four bosons. Hence there are 
some terms which are not created so that they do not 
need to be neglected in the truncation scheme. But this 
effect is obviously fairly small for the system under study. 

From these observations we will restrict our compu- 
tations for the QO to the schemes which are based on 
stationarity, i.e., we do not show the direct Hartree-Fock 
calculation. We performed such calculations and the re- 
sults were unsatisfactory. The flow equations did not 
even converge to a diagonal Hamiltonian fixed point for 
A > 4. This illustrates again that the starting point 
should not be too far away from the exact result^ 

Furthermore, we show results only for the MKU- 
generator. Its results are representative for similar calcu- 
lations based on Wegner's generator. They are shown in 
Fig. 0] Moreover, the MKU-generator is easier to imple- 
ment (compare Eqs. Ijl7|) and Ijl8(l above). This fact will 
be of importance for calculations comprising more terms 
(see below). 

The accuracy of the results is very impressive in view 
of the simplicity of the numerical effort. As for the PQO 
the variation of the renormalized, i.e., final ground state 
energy yields the best results if the focus includes the 
two-particle energies whose accuracy is given by r-i. 

Finally, a remark about another possible choice of the 
generator is in order since its generalisations have been 
used in more complex physical situations, e.g. in Refs. 
fin! Using the W-VAR-^ 1 q scheme, the quadratic 
off-diagonal coefficient 32,0 vanishes at I = 0. But it 
will become finite in the course of the flow. It is pos- 
sible to modify Wegner's generator to keep this coeffi- 
cient zero during the flow. To this end, the 772,0 co- 
efficient of the generator is changed. It is chosen to 
be r7 2 ,o = -7252,134,0(50,1 + 33o,2)/(ffo,i + 90,2 + 634,0) 
instead of 2(go,i + 30,2)32,0- This choice has two ma- 
jor advantages. First, the coefficient 32,0 does not need 
to be traced. Second, we find that it converges faster 
to its diagonal Hamiltonian. For the PQO, the accu- 
racies are found to be 10 2 r = 0.00569, 10 2 ri = 1.53 
and 10 2 7*2 = —13.6. Thus we obtain the best value of 
the ground state energy with this scheme but also the 
worst value of the single-particle energy. The value of 
the interaction is as inaccurate as it is for the W-VAR- 
3g Q scheme. Furthermore, we have lost the possibility to 
vary the initial value of the parameter a since it is fixed 
to unity by the condition (72.0 = 0- An adapted choice of 
a will be important in dealing with more terms so that 




2 4 6 8 10 

A 



MKU-VAR-flS' MKU-VAR-sS' 

MKU-VAR-gg' MKU- VAR-g^' 



2 4 6 8 10 

A 

FIG. 4: Relative errors ro, ri and 7*2 of the first three diag- 
onal coefficients g^j , g^± , g^ ( see text) for four variation 
schemes based on the MKU-generator. 



we consider it a caveat that a is already fixed. Summa- 
rizing, this variant yields very reasonable results, but it 
is not the optimum choice. 



2. Spectral Weights 

Apart from the spectrum of the oscillator, one is also 
interested in the computation of the Green function and 
of the spectral weights. The retarded Green function is 
defined as 

G(t) = -i(0\[X(t),X(0)}\0)e(t), (20) 
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where X(t) = exp(iHt)X exp(— iHt), |0) is the ground 
state. The Fourier transform of the Green function reads 



G(w) 



+ O0 



dtG(t) exp(iwi) 

) 



1 



(21) 
(22) 



uj - (E n ~ E ) + 
1 

uj + (E n - E ) + i0+ 



We have denoted the eigenstates of the Hamiltonian by 
\n) with eigen energies E n , i.e. H\n) = E n \n). We 
fill focus on the computation of the spectral weights 
J n := |(rt|X|0)| 2 . They can be found trivially after the 
transformation since then the eigen states are the states 
with a given number of bosons 



|0) = [/(oo)|0) 
|n) = J7(oo)|n) 

= C/(oc)V>>o 



(23a) 
(23b) 

(23c) 



So we need the transformed observable X(l) = 
W(l)XU(l) which satisfies the same flow equation as the 
Hamiltonian 



d l X(l) = [r,(l),X(l)] 



(24) 



with initial condition X(0) — (<r/\/2) (a^+a). Here again, 
an infinite hierarchy of terms appears when X(l) is com- 
muted with r){l). We will restrict the calculation to the 
lowest non-trivial truncation scheme 

X(l) = aifi(l) (a 1 * + a) + a 1A (l) (a) 2 a + a f a 2 ^ 

+a 3 ,o(0(a t3 + « 3 ) ■ ( 25 ) 
The flow equations of the a coefficients are found to be 

diaifl = -2[7? 2 ,o(Q!i,o + ais + 3a 3 , ) (26a) 
+3?72.ia3,o + 12»y4,oa3,o] 

diai.i = -2[?72,o(2ai,i + 3a 3i o)] (26b) 
-3[??2,i(ai,o + 2(21,1 + 6a 3 , ) + 12774, a 3 ,o] 

dia 3 . Q = -2ry 2 ,oo:i,i + ?72,iai,o (26c) 
-4??4,o(ai,o + 301,1) ■ 

Note that this set of equations is linear in the a coef- 
ficients since the coefficients in the generator are deter- 
mined from the flow of the Hamiltonian. So for the com- 
putation of the observables the g coefficients are given 
and only a linear set of differential equations, though 
with non-constant coefficients, has to be solved. 

The first three J coefficients can be obtained from (|26[) : 
they read J\ = a\ o(°°)? J2 — due to the identical par- 
ity of the |0) and |2) states, and J 3 = 6a 3 (oo). The 
parity of the states in real-space representation is ex- 
pressed in second quantization by the fact that the boson 



number is changed only by multiples of two. In Fig. |SJ 
the results obtained in this way, for the MKU generator, 
are compared to the numerically exact ones. The rela- 



tive errors S1/3 = J]y 3 — /J1/3 arc plotted. For 



7-(° x ) 
1/3 



10 3 Si 





1 1 


1 1 




















1/ 


MKU-VAR-gS' - - - 


-■ MKU-VAR-g^' - 


1/ 

1/ 


MKU-VAR-gg' 

1 1 


- MKU-VAR-gW 

1 1 



10 



10s 3 





1 1 

MKU-VAR-flS' --- 


1 1 

-- MKU-VAR-45,' 






^MKU-VAR-gg' 


■■- MKUAAR-gg' — 






















1 1 


1 1 







2 4 


6 8 


: 



10 



FIG. 5: Relative errors si and S3 of the CUT calculation for 
the spectral weights Ji and J3 (see text) for four variation 
schemes based on the MKU-generator. 

A = 10, the error in J\ is of the order of a few tenth of 
a percent. The error on J 3 is about 10%. Note that the 
weight Ji (j} ex) (A = 10) = 1.297 • 10" 1 ) dominates and 

that the weight of J 3 (J^° x) (A = 10) = 3.469 • 10" 4 ) is 
smaller by more than two orders of magnitude so that 
the errors on the absolute scale are equally small. We 
attribute the relatively large deviations in J 3 to the low 
order of the truncation scheme used for the observables 
and for the Hamiltonian where no three-particle energy 
was included. 

For the PQO, the relative errors of the computation of 
the spectral weights are given in Table ITTA for MKU and 
W generators. The exact values for the spectral weights 
are j[ cx) = 0.521 and J^ cx) = 0.00152, showing that even 
for the PQO the lowest excited state has most of the 
spectral weight. 

It is very interesting to compare the various varia- 
tional schemes. The result from the energies is con- 
firmed. The by far best result is obtained from the vari- 
ation of the renormalized ground state energy. This fact 
strongly suggests to adopt this scheme also in more com- 
plex situations. The variation of the bare ground state 
energy, equivalent to the CUT treatment around the self- 
consistent mean-field calculation, yields also satisfactory 
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10 a si 


10s 3 


W-VAR-3^ 


2.77 


-1.99 


W-VAR- fl ^ 


2.02 


-1.00 


W-VAR-3^ 


3.58 


-1.15 


W-VAR- fl ^ 


3.96 


-1.24 


MKU- VAR- gjffl 


3.71 


-2.94 


MKU-VAR-g^/ 


2.12 


-1.23 


MKU- VAR- 


3.17 


-1.46 


MKU-VAR- fl ™ 


4.06 


-1.80 



TABLE III: Relative errors Si and S3 of the CUT calcula- 
tion for the spectral weights Ji and J3 of the PQO for four 
variation schemes based on the W and MKU generators. 



results. It is easier to realize which may become essential 
for extended systems. Finally let us note that both gen- 
erators give results which are similar, as was already the 
case for the spectrum. Wegner's generator works a little 
better for the spectral weights than the MKU generator, 
whereas it was the inverse for the spectrum. 



IV. EXTENDED TRUNCATION SCHEMES 

In this section, we only consider the PQO of Eq. 10 to 
keep the presentation concise. The question investigated 
is how the truncation scheme used so far can be extended 
to improve the results for the spectrum in a systemati- 
cally controlled fashion. For simplicity, we stick to the 
MKU-generator. 



A. Structure of the Flow Equations 

Using the MKU-generator, the Hamiltonian retains its 
band structure during the flow so that no terms changing 
the boson number by more than four arise. Hence the 
generalisation of Eq. (|12fl reads 



H 



g ,iM l + g2,ia {2 M l + gi^Mia 



(27) 



4.4 4 
+34, i a ^ + g^iMid 



where summation over the repeated index i from to 00 
is understood. The same convention is used in the sequel 
unless upper bounds for the indices are explicitly dis- 
cussed. The operator M$ — a^a appearing in Eq. ljT7|) 
is the normal-ordered form of the i th power of the num- 
ber operator Q = a^a. In Wegner's scheme the flowing 
Hamiltonian would contain terms creating or destroying 
any even number of the excitations. The MKU generator 
is given by 



gi,iMia + g iti 



a) A Mi - g^Mia 



digo,k = - 


0O.2.2 

Pk,i,j 92,i92,j " 


a 0,4,4 

^Pk,i]j 9i,i9u 


di92,k = - 


a 2,0,2 

"k,i,j 90,i92,j - 


a 2,2,4 
1- P k,i,j 92,i9ij 


digi.k = - 


a 4.0.4 

P k,i,j 9o,i9i,j 





The flow equations are found by commuting 77 MKU with 
H . The derivation of these equations is included in Ap- 
pendix [5] where also the relations needed for Wegner's 
scheme are given. The flow equations have the general 
form 



(29a) 
(29b) 
(29c) 



Let us for example consider the g$s which appear in the 
part of the Hamiltonian conserving the number of exci- 
tations. They are renormalized only by the commutation 
of one operator creating n = 2 (4, resp.) excitations and 
of one operator destroying n = 2 (4, resp.) excitations. 



B. Truncation Schemes 

In the numerical integration of the differential equa- 
tions, we can keep only a finite number of couplings. We 
denote by mo, m,2 and 7714 the number of 170 , 52 and 34 
couplings kept. For the go couplings, this means that 
we keep the couplings 170,0, ■ ■ ■ , 3o,m -i- We will call the 
truncation scheme defined by the three numbers mo, TO2 
and TO4 the "(too,to,2,TO4) scheme". For example, the 
truncation scheme used in Sect. IIII C II is the (3,2,1) 
scheme. 

One question arising is how to choose the trunca- 
tion scheme. This is a very difficult question which to 
our knowledge has no general answer. Starting from a 
Hamiltonian of the form (mo, 7712, TO4) and commuting it 
once with the corresponding MKU-generator, new terms 
are generated which belong to the (max[2m,2, 2(to4 + 
1)], max [mo + — 2, TO2 + rrii], mo + TO4 — 2) scheme. 
Starting from the initial (3, 2, 1) scheme, for instance, the 
following schemes are generated 



(3,2,1) 



+ (4,3,2) 
(10,9,8) 



(6,5,4) 
(18,17,16) 



(30) 



(28) 



These schemes are all of the type (N,N - 1,N -2). So 
we decide to consider the schemes of this family. Addi- 
tionally, the terms of this family have the nice property 
that in each class changing the number of excitations Q 
by 0, 2 or 4 the maximum number of creation and an- 
nihilation operators is equal to 2(N — 1) = 2 (too — 1) = 
2 + 2(m 2 - 1) = 4 + 2(to 4 - 1). 



C. Results 

We have checked that including more terms in the 
Hamiltonian, i.e., increasing N in the (N, N — 1, N — 2) 
scheme, results in a continuously improved accuracy for 
the low-lying energy levels. In Fig. we illustrate re- 
sults for the ground-state energy on passing from the 
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scheme (4,3,2) to the scheme (8,7,6). For both trun- 
cation schemes, the ground-state energy is plotted as a 
function of the parameter a (see Sec. IHIjl . Clearly, a 
plateau behaviour develops around a value of a smaller 
than 1. Note in particular that the magnitude of the scale 
on which the ground state energy changes is shrinking 
strikingly from the (4, 3, 2) to the (8, 7, 6) scheme. This 
finding agrees perfectly with our previous argument on 
the independence of an exact calculation on the start- 
ing point, see Sect. 1111 Bl Hence an improved calcula- 
tion should approach this independence gradually. This 
means that the dependence of the quantities computed 
decreases as the accuracy of the calculational schems is 
enhanced. This is exactly what we can observe. The fact 
that values of a below 1 are preferred is plausible. The 
X 4 part of the potential squeezes the eigen states more 
than does the harmonic potential. 




0.5 1.0 1..5 



a 



1 1 1 1 1 


(8,7,6) 


So°V) G [0.36760693,0.36760701] 


i 


i i i i 



0.4 0.8 1.2 



a 

FIG. 6: Ground state energy of the PQO as a function of 
a for the two schemes (4, 3, 2) and (8, 7, 6). The dashed line 
represents the exact ground state energy. 

Fig. shows also another feature. The results deterio- 
rate if a becomes too big. In fact, the flows even diverge 
above some critical value of a. For the (4, 3, 2) scheme 
this value is er c ~ 1.52 and for the (8,7,6) scheme it 
is fi c ~ 1.21. These divergences originate in the error 
induced by the truncation of the flow equations. The 
values of a c decrease on increasing N. For the values we 
investigated (up to 25) the values of a c remained finite. 

These results suggest that there is no uniform conver- 
gence as the limit N — > oo is taken. For any practical 
calculation, however, this is no serious limitation since 
the values of N cannot be chosen very large. In particu- 



lar in spatially extended systems, N — 3 or 4 will be the 
maximum which can be treated. 



V. SUMMARY AND CONCLUSION 

In this article, we presented a detailed investigation of 
the quantum quartic oscillator by perturbative and non- 
perturbative, renormalizing treatments. We applied con- 
tinuous unitary transformations to determine the spec- 
trum as well as the spectral weights of the quartic os- 
cillator. This model served as the simplest interacting 
quantum model. So the method could be illustrated in 
a simple and transparent setting. The calculations pre- 
sented may guide future investigations of spatially ex- 
tended systems like, for instance, interacting spin waves 
in quantum antiferromagnets. 

We demonstrated that it is rewarding to ally the CUTs 
with a variational principle. Thereby, an optimum refer- 
ence state can be chosen with respect to which the neces- 
sary truncations introduce only very small errors. It is to 
be expected that this is also true for other renormalizing 
schemes. 

In the CUT treatment presented tracing only six cou- 
plings achieved a remarkable accuracy. Without a rea- 
sonably chosen reference state, however, the CUT pro- 
cedure does not succeed. Information on the spectral 
weights complements the information needed for the 
Green functions. Only little additional effort is needed 
to obtain the required coefficients. 

We compared two generator schemes (Wegner (W), 
Eq. J2J and Mielke-Knetter-Uhrig (MKU), Eq. Q). The 
W generators always lead to a fixed point while the MKU 
generators may fail to converge. But the W generators 
are stopped by degeneracies so that no significant sim- 
plification is reached. The MKU generators are easier 
to write down and to implement than the W genera- 
tors. In their performance they are similar for the quar- 
tic oscillator. For the energies the MKU is slightly better 
than the W generator. For the spectral weights it is the 
inverse. The experience with more complex, extended 
systemsi&iS indicates that the MKU generator is very 
powerful if the energy eigenvalues are sorted in ascend- 
ing order of the particle number Q. 

Finally, we discussed how the truncation scheme for 
the continuous unitary transformation can be extended 
to improve the accuracy of the approach in a systemat- 
ically controlled way. Our results show that this is in- 
deed possible. Although no uniform convergence can be 
expected the description of the low-lying states can be 
gradually improved by keeping more terms in the flow. 

Ideally, our presentation will serve as seed for future 
investigations of more complex systems by the approach 
discussed here. 
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couting arguments, and one finds {Mi — a) % a) 



a , a' 



ip-n 
i 1 



£ (c\C c n C;) M„_ C) (Bl) 



= n\/\p\(n — £>)!] being the usual binomial coefficient 
for p ^ n, and being otherwise. The sum in 1|B1|) 
is performed over the number c of contracted a and 
operators. 



APPENDIX A: PROOF OF CONVERGENCE OF 
WEGNER'S FLOW FOR INFINITE MATRICES 



2. Derivation of the Flow Equations 



The proof of convergence given in Ref. [lj is only valid 
for finite matrices since it makes use of the invariance 
of the trace of H 2 (l) in the flow. But the trace is not 
generally defined for infinite systems. So the proof as 
given in Ref. Q works for the two- level system l|5a|) , but 
it does not work for the bosonic oscillator (|5b|) whose 
eigenvalues grow like n. 

But the proof can be slightly modified so that it works 
for infinite matrices as well. The idea is inspired from 
statistical mechanics. A regularization is needed which 
makes the higher lying eigen values contribute less. To 
this end, one can study exp(—(3H(l)) instead of H for an 
arbitrary but fixed j3 > 0. Since exp(—/3H) and H com- 
mute there is a certain basis in which both are diagonal. 
So the aim is to diagonalize exp(— (3H (I)) which guaran- 
tees the diagonality of H up to possible degeneracies. 

The generator is modified 



1 



m (l) = [exp(-/3ff(Z))| d ,exp(-/3ff(0)| nd ] 



(Al) 



Then Wegner's proof can be used with H replaced by 
exp(— /3ff). The division by /3 2 in r\ only changes the scale 
of I. But it is useful since it ensures lim 1 3_ > o ^ = %r. 
Thus one recovers in this "high-temperature" limit the 
generator proposed by Wegner. Hence the flow equation 
diexp(—(3H(l)) = [r)p(l),exp(—{3H(l))] can be reduced 
to the usual one diH(l) = [r] w (l), H(l)] because it does 
not matter whether the unitary transformation is applied 
to the whole function or to its argument. 

So the use of (5 > serves only to induce convergence. 
For (3 — > all states of the Hilbert space matter as before. 



APPENDIX B: FLOW EQUATIONS 

1. Basic Relation for Commutators of Bosonic 
Operators 



Here the evaluation of the commutator of the generator 
with the Hamiltonian is addressed. Note that for Weg- 
ner's generator the first task is to compute the generator, 
given by the commutator of two parts of the Hamilto- 
nian. As the computation of r] w is qualitatively not dif- 
ferent from the computation of [rj. H] we will not discuss 
it here but focus on the MKU-generator. We will how- 
ever use notations which are more general than the ones 
used in the main text, and that are easily transferable to 
Wegner's generator. 

The flowing Hamiltonian has the form 



E 

k=2A 



H = H 



with H~ = (H+y and 



H a = Yl 9o,i M i and H k 

iGN 



ieti 



(B2) 



(B3) 



For Wegner's generator, the sum over k in i|B2D would 
extend over all non-negative even integers. The MKU 
generator simply reads 



7?MKU 



2^ W 

fe=2.4 



E 

k=2A 



(Ht-H k 



from which one deduces the flow equations 

diH = 2([H+,H^] + [H+,H^]) 
d t H+ = [H+,H ]+2[H+,H^] 
diHt = [H+,Ho\. 



(B4) 



(B5a) 
(B5b) 
(B5c) 



As an example, let us now outline the computation of 
the commutator [H^ - , H~] for n ^ p. It is equal to 



2^ 9 P ,i'9 



ip-n 
= a' 



E 9p,i'9n,i" E A id" - Ml ' 



(B6a) 

" : n c (B6b) 



The computation of the commutator of the genera- 
tor with the Hamiltonian or the observables basically 



involves commutators 



It suffices to consider 



n p since if n ^ p, one can use the relation [a™, aJ P ] 



a , a 1 



Such a commutator is computed using basic 



,t*- 



"E^ 

i 

^2^-i'%i'-n+c,c9p,i'9n,i-i'-n+c , (B6c) 
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where we defined 



(B7) 



In (IB6c|l . the i runs over N, but the sums over i' and c 
are performed on a finite set of positive numbers, because 
C% = if p > n. Note that by setting n = in HB6cf) . 
one obtains [JT+,i?o] 

Collecting everything gives 



+-^i'\i-i'-4 + c,c9i,i'94,,i-i'-i+c (B8a) 



+ 2 - 4 ^-i'-2+ C ,c5'4, J '52 :l -i'-2+ C J (B8b) 



The flow of the observable can be computed in the 
same way. One has 



x = E 



(B9) 



with = and 



The flow equations are found to be 



diX+ 
d t X+ 



[H+,Xi] - [H£,X+] 
+ [H+,Xz] - [HZ,X+] 
[H 2 ,X X ] - [H^,X h ] 
+ [H+,X-] [H-,X+] 

+ [ H t' X k-A\ ~ l H i ' X fc+J 



(BIO) 



(Blla) 
(Bllb) 
(Bile) 



These equations can be written in an explicit form, as we 
did in transforming l|B5(l into JEHJ. But we refrain from 
doing so, since they would be more complicated than 
(|B8|) . All relations needed to derive them is contained in 
this Appendix. 



/c£2N+l 
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diagonal part of the Hamiltonian and the X 2 term as the 
non-diagonal part. Writing H(l) = [a(l)P 2 + f3(l)X 2 ]/2, 
it is found that 77W = — 2iaf3(XP + PX) so that the flow 
equations dm = 2a 2 [3 and dia — —2a 2 (3 result. Thus a/3 
is a constant, and as (3 — > 0, a — > 00. Of course, this was to 
be expected since the spectrum of an harmonic oscillator 
is discrete so that it cannot be mapped one-to-one to the 
continuous spectrum of P 2 . 



